Low-energy resonances and bound states of aligned bosonic and fermionic dipoles 
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The low-energy scattering properties of two aligned identical bosonic and identical fermionic dipoles are 
analyzed. Generalized scattering lengths are determined as functions of the dipole moment and the scattering 
energy. Near resonance, where a new bound state is being pulled in, all non-vanishing generalized scatter- 
ing lengths diverge, with the aoo ar, d «n scattering lengths being dominant for identical bosons and identical 
fermions, respectively, near both broad and narrow resonances. Implications for the energy spectrum and the 
eigenfunctions of trapped two-dipole systems and for pseudo-potential treatments are discussed. 
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Currently, the creation of ultracold heteronuclear ground 
state molecules poses one of the major experimental chal- 
lenges in the field of ultracold physics yj]. The trapping of 
ultracold ground state molecules with large phase space den- 
sity promises to allow an exciting array of novel research lines 
to be studied. Although the largest phase space density of ul- 
tracold ground state molecules achieved to date is still fairly 
small, a number of promising cooling schemes have been 
demonstrated ]2|] . Thus, it is expected that degenerate molec- 
ular gases with large electric dipole moment will be created in 
the laboratory in the near future. Polar molecules are a can- 
didate for qubits in quantum computing [3] and may be used 
in high precision measurements that aim at placing yet stricter 
limits on the electric dipole moment of the electron jj]. Fur- 
thermore, dipolar gases are predicted to show roton-like fea- 
tures |H1] and to exhibit rich stability diagrams whose details 
depend on the trapping geometry The stability of dipolar 
atomic Cr condensates has recently been investigated exper- 
imentally. To enhance the anisotropic effects, which are due 
to Cr's magnetic dipole moment, the s-wave scattering length 
was tuned to zero by applying an external field in the vicinity 
of a Fano-Feshbach resonance (7|]. 

To create and then utilize ultracold molecules, it is manda- 
tory to develop a detailed understanding of the scattering 
properties of two interacting dipoles in free space and in a 
trap. Unlike the interaction between s-wave alkali atoms, 
the interaction between two dipoles is long-range and angle- 
dependent. A two-dipole system can, e.g., be realized experi- 
mentally by loading an optical lattice with either two or zero 
dipoles per site. If the optical lattice is sufficiently deep and if 
the interaction between nearest and next to nearest neighbors 
are absent or negligible, then each optical lattice site can be 
treated as an independent approximately harmonic trap. 

This paper determines the scattering properties of two 
aligned dipoles, either identical bosons or identical fermions, 
as functions of the dipole moment and the scattering energy. 
In general, the dipoles can either be magnetic or electric. For 
concreteness, we restrict our discussion in the following to 
the scattering between molecular electric dipoles. Sequences 
of scattering resonances, which can be classified as "broad" 
and "narrow", are found. For identical bosons, these reso- 
nances have previously been termed potential and shape res- 
onances, respectively, and have been interpreted within the 
framework of adiabatic potential curves [8]. The resonance 
positions are correlated with the appearance of bound states in 



free space and "diving" states in the energy spectrum of two 
aligned dipoles under external confinement. The nature of the 
broad and narrow resonances is further elucidated by analyz- 
ing the bound state wavefunctions. In addition, we show that 
the eigenequation of two aligned dipoles in a harmonic trap in- 
teracting through an anisotropic zero-range pseudo-potential 
reproduces much of the positive energy spectrum, but exhibits 
some peculiar unphysical behavior for small energies. The 
origin of this unphysical behavior is pointed out and a simple 
procedure that eliminates it is presented. 

Neglecting hyperfine interactions and treating each dipole 
as a point particle, the interaction potential between two 
dipoles aligned along the z-axis is for large interparticle dis- 
tances r given by Va(jr), V&i(r) = d 2 (l — 3cos 2 9)/r 3 . Here, 
d denotes the dipole moment and the angle between the rel- 
ative distance vector f and the z-axis. We model the short- 
range interaction V sr (r) between the dipoles by a simplistic 
hardwall potential, V sr (r) = °° for r < r c and for r > r c , so 
that the full model potential is given by V m (r) — V ST (r) for 
r < r c and Va{f) for r > r c . The boundary condition im- 
posed by r c can be thought of as introducing a short-range 
/iT-matrix, which is modified by the long-range dipole poten- 
tial @1 . The characteristic length scale of V sr (r) is given by 
the hardcore radius r c and that of Vdd(r) by the dipole length 
D*,D*= /jd 2 /h 2 , where /j denotes the reduced mass. The cor- 
responding natural energy scales are given by E,- c and Ed„ , re- 
spectively \E rc = h 2 /(fir 2 ) and Ejj t = fr 2 /{jiD 2 )]. A straight- 
forward scaling of the relative Schrodinger equation shows 
that D t and r c are not independent but that the properties of 
the system depend only on the ratio D*/r c lloll . This ratio can 
be tuned experimentally by varying D t through the applica- 
tion of an electric field |9[] . 

To obtain the ^-matrix elements K^'™ 1 ' , where I and I' de- 
note the orbital angular momentum quantum number of the 
incoming and outgoing partial waves, respectively, and »?/ and 
my the corresponding projection quantum numbers, we solve 
the relative Schrodinger equation for V m (r) for a fixed scat- 
tering energy E sc numerically. The azimuthal symmetry con- 
serves the projection quantum number, and throughout we re- 
strict our analysis to m/ = 0. The radial Schrodinger equation 
is propagated using the Johnson algorithm with adaptive step 

size fiUl . The ^-matrix elements K l ,Q(k) = tan8//(fc) are 
found by matching the log-derivative to the free-space solu- 
tions at sufficiently large r. Since the long-range part of V m (r) 
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Figure 1: Scaled scattering lengths (a) aoo Ac > (b) a 2o/ r c and (c) 
a 22 l r c as a function of the scaled dipole length D*/r c for two identi- 
cal bosons interacting through V m for three different scattering en- 
ergies: E sc = 9.36 x ltT 8 ^ (solid line), E sc = 9.36 x 10~ 6 E rc 
(dashed line) and E sc = 9.36 x l0~ s E rc (dotted line). Crosses and 
squares in the top panel indicate the positions of the broad and nar- 
row resonances, respectively, predicted by analyzing the WKB phase 
of the adiabatic potential curves (see text). 



is proportional to the spherical harmonic Y2o(Q,§), the phase 
shifts $ijr(k) are only non-zero if \l — 1'\ < 2. 

Figures Q] and [2] show the generalized scattering lengths 
ai ii for two identical bosons and two identical fermions, re- 
spectively, for three different scattering energies E sc as a 
function of the dipole length D t . The scattering lengths 
aijr(k), a/ [/(k) = —Kiji(k)/k (k denotes the wavevector, k = 

are defined so that the a/ /' (k) approach a con- 
stant as k — > lfl2l[T3ll . The largest D*/r c value considered in 
Fig. Q]is 40. If we choose r c « 10ao, then the largest dipole 
length considered in Figs. Q] and [2] is D™ ax w 400ao, imply- 
ing a minimum dipole energy £™ n of 1.27 x 10~ 4 K. For the 
polar molecule OH, this corresponds to a maximum dipole 
moment of 0.404 Debye, a value that should be attainable ex- 
perimentally. The scattering energies in Figs. Q] and [2] range 
from 9.36 x 10~ 8 £V C to 9.36 x 10~ 5 £V f , or, using as before 
r c = lOflo, from 1.91 x 10~ 8 K to 1.91 x 10~ 5 K. Thus, the 
largest E sc /Eo t value considered in Figs. Q]and[2]is 0.15. This 
places the present study in the regime where the minimum 
value of the cross section has been predicted to behave univer- 
sally Hill , but where the parameters of the two-body potential 
and the s-wave scattering length it results in, especially near 
resonance, are important 11511 . 

FigureQJa) shows the scattering length aoo as a function of 
for two identical bosons interacting through V m (r). Five 
broad and two narrow resonances (located at w 23r c and 
37r c ) are clearly visible. FiguresQlb) andQJc) show the gener- 
alized scattering lengths 020 and 022, respectively. In the Born 
approximation (BA) for V^r), both 020 and 022 vary linearly 
with lfi"3l 12111 . 020 and «22 obtained from the full coupled 



Figure 2: Scaled scattering lengths (a) a\\/r c and (b) a^i/r c as a 
function of the scaled dipole length D*/r c for two identical fermions 
interacting through V m for the same three scattering energies as in 
Fig.frj Crosses and squares in the top panel indicate the positions of 
the broad and narrow resonances, respectively, predicted by analyz- 
ing the WKB phase of the adiabatic potential curves (see text). 



channel calculation show deviations from the BA for certain 
D t values. The positions of the "spikes" coincide with the res- 
onance positions of aoo- Notably, the widths of the spikes de- 
crease with increasing l + V , The top panel of Fig.Q]shows the 
resonance positions as predicted by the WKB phase accumu- 
lated in different adiabatic potential curves |8|]. The crosses, 
obtained by analyzing the WKB phase of the lowest adiabatic 
potential curve, predict the positions of the broad resonances 
very accurately. The squares, obtained by summing the WKB 
phases of all other adiabatic potential curves, predict the num- 
ber of narrow resonances semi-quantitatively but do not pre- 
dict their positions accurately . 

Figures |2ja) and |2jb) show the generalized scattering 
lengths an and 031 for two aligned identical fermions inter- 
acting through V m (r) as a function of D*. Away from res- 
onance, an and 03 1 vary approximately linearly with D t . 
The spikes in Fig. [2] are interpreted as resonances, which 
we term, as in the boson case, broad and narrow B22I1 . Fig- 
ure [2] shows five broad and two narrow resonances (located 
at D* w 23.5r c and 37.5r c ). A key difference between dipole 
scattering of identical bosons and identical fermions is that the 
lowest non-vanishing scattering length for bosons (i.e., aoo) 
cannot be approximated by applying the BA to V^r) (the 
BA for Vdd(r) gives aoo = 0) while the lowest non-vanishing 
scattering length for fermions (i.e., an) can be, away from 
resonance, approximated by the BA for Vddt/) (the BA for 
Vdd(r) gives an = — 2D*/5) 11131, |2 ill . The crosses and squares 
shown in the top panel of Fig. [2] indicate the positions of the 
broad and narrow resonances, respectively, as predicted from 
the WKB phase of the lowest adiabatic potential curve and of 
all other adiabatic potential curves. For identical fermions, the 
WKB prediction for the positions of the broad resonances is 
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Figure 3: Eigenenergies for two identical bosons in the vicinity 
(a) a broad resonance and (b) a narrow resonance, and for two iden- 
tical fermions in the vicinity of (c) a broad resonance and (d) a nar- 
row resonance. Solid lines and circles show the bound state ener- 
gies for two dipoles in free space interacting through V m and V ps , 
respectively. Dashed lines and crosses show the bound state ener- 
gies for two dipoles under external harmonic confinement interacting 
through V m and Vp S , respectively. The oscillator length flho is given 
by y/h/JJIa), and the hard core radius r c of V m is 0.00306 ciho- 



Figure 4: Eigenfunctions for two identical bosons in the vicinity of 
(a) a broad resonance (D* = 0.1061aho) and (b) a narrow resonance 
CD* = 0.07147ahoX and for two identical fermions in the vicinity of 
(c) a broad resonance (D* = 0.07976aho) and (d) a narrow resonance 
(D* = 0.07167aho) m a spherical harmonic trap interacting via V m for 
9 = 0° (solid), 9 18.2° (dashed), 9 36.4° (dotted) and 9 w 54.6° 
(dash-dotted line). In all panels, r c equals 0.00306 a^ D . Note the log 
scale for r and the different ranges of the y-axis. 



less accurate than that for identical bosons. 

Figures Q] and [2] show that the widths of broad and narrow 
resonances increase with increasing E sc for fixed D*/r c and 
with increasing D*/r c for fixed E sc . Putting this together, we 
find that the resonance widths of both broad and narrow reso- 
nances increase with increasing E sc /Ed^. 

To better understand the resonance structure in Figs. Q] 
and|2] we determine the bound state energies of the two inter- 
acting dipoles in free space. The Schrodinger equation for the 
relative coordinate is solved using two-dimensional B-splines. 
The two-dipole system supports a new bound state at those 
D*/r c values where the scattering lengths aoo and an for two 
identical bosons and fermions, respectively, diverge. Solid 
lines in Figs. [3ja) and |3jb) show the bound state energy for 
two identical bosons in the vicinity of a broad and a narrow 
resonance, respectively, while solid lines in Figs. 0c) and[2d) 
show the bound state energy for two identical fermions in the 
vicinity of a broad and a narrow resonance, respectively. 

The two-body energy E/, of weakly-bound s-wave interact- 
ing systems is well described by the s-wave scattering length 
aoo, Eb = — h 2 / (2fj{aQo(E)] 2 ) . To test if this simple pseudo- 
potential expression holds for dipolar systems, we analytically 
continue the scattering lengths for V m (r) to negative ener- 
gies. We obtain stable aoo(E) and an(E) for negative scat- 
tering energies by matching the coupled channel solutions to 
the free-space solutions at relatively small r values (r max w 
The bound state energies for two identical bosons 
in free space, determined self-consistently iflrSl Il7tl from 
Eh = — fir / (2fj{aQo(Eb)] 2 ), are shown by circles in Figs.[3ja) 
and Ob). Similarly, we determine the bound state energies 
for two identical fermions in free space by solving the equa- 
tion Eb = —h 2 /(2fj{au(Ei,)] 2 ) 0181H9I1 self-consistently [cir- 
cles in Figs. |3lc) and[3d)]. Somewhat surprisingly, the bound 
state energies in the vicinity of both broad and narrow reso- 



nances are very well described by a single-channel expression 
for identical bosons and identical fermions (see below for fur- 
ther discussion). 

In addition to the free-space system, we consider the 
trapped system. Dashed lines in Fig. [3] show the energies 
for two dipoles interacting through V m (r) under external har- 
monic confinement Vt ra p, Vtrap = /jcoV 2 /2. Near resonance, 
the lowest state with positive energy changes rapidly and turns 
into a negative energy state with molecular-like character. The 
energy of this "diving" state is slightly higher than the energy 
of the free-space system (the trap pushes the energy up). 

Figures Ufa) and |4|b) show the scaled eigenfunctions 
r\\t(r, 0) for two identical bosons with E rs — 0.88fico as a func- 
tion of r for different near a broad and a narrow resonance, 
respectively. Similarly, Figs. Hfc) and |4jd) show the scaled 
eigenfunctions for two identical fermions with E w — 0.88/ico 
near a broad and a narrow resonance, respectively. In all pan- 
els, the wavefunction cut for = 0° has the largest amplitude, 
reflecting the fact that the dipole-dipole potential is most at- 
tractive for = 0°. Interestingly, the nodal structure of the 
wavefunction for two identical bosons near a broad resonance 
[Fig.HJa)] has a similar structure to that of the wavefunction 
of two identical fermions near a broad resonance [Fig.|4fc)]: 
Both nodal surfaces show approximately spherical symmetry. 
On the other hand, the nodal structure of the wavefunction for 
two identical bosons near a narrow resonance [Fig. |4jb)] has a 
similar structure to that of the wavefunction for two identical 
fermions near a narrow resonance [Fig.|4fd)]: The nodal sur- 
faces depend on both r and 0. To quantify the higher partial 
wave contributions, we project the wave functions shown in 
Fig- H] onto spherical harmonics. The s-wave contribution of 
the boson states near the broad and narrow resonances is about 
95%, while the p-wave contribution of the fermion states near 
the broad and narrow resonances is about 95 and 80%, re- 
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spectively. The gas-like states near resonance, in contrast, are 
dominated by a single partial wave (for bosons, e.g., the s- 
wave contribution of the energetically lowest-lying gas-like 
state is about 99% while the of-wave contribution of the ener- 
getically next higher-lying state is about 99%). In the future, it 
will be interesting to investigate how the higher partial wave 
contributions of the weakly-bound anisotropic molecules af- 
fect the scattering properties between two such composite 
particles and, more generally, the BEC-BCS crossover-type 
physics. 

Lastly, we show that the entire energy spectrum of two 
aligned dipoles under external harmonic confinement interact- 
ing through V m (r) can be reproduced by a zero-range pseudo- 
potential framework. Since the anisotropy of the dipole-dipole 
interaction leads to a coupling of different partial waves, 
the pseudo-potential Vp S (r) contains an infinite number of 
terms HcJ, V ps (r) = Y!j°ji = ogi,i'{k)®ij>{r), where the coupling 

strength gw(k) is proportional to — tanS/ p {k) / k l+l ' +l and 
©;//(r) denotes an operator 1I20I1 . Assuming the dip vanish 
for \l — 1'\ > 2, as is the case for two interacting dipoles, the 
eigenequation for two particles under spherically symmetric 
external harmonic confinement can be elegantly written in 
terms of a continued fraction [21]. To obtain the eigenen- 
ergies for two aligned dipoles under external harmonic con- 
finement interacting through V ps (r), we solve the eigenequa- 
tion self-consistently, using the energy-dependent a; // (k) ob- 
tained for V m (r) as input parameters. The resulting ener- 
gies, shown by crosses in Fig. [3] agree well with those ob- 
tained for V m (r) (dashed lines). However, for small \E\ the 
eigenequation for the pseudo-potential results in an unphysi- 
cal eigenenergy [not shown in Figs. Eta)-(d)]. For two iden- 
tical bosons, e.g., the eigenequation for V vs (r) permits a so- 
lution with E w 0.05to, which is absent in the eigenspec- 
trum of two identical bosons under external harmonic con- 
finement interacting through V m (r). Importantly, if we restrict 
the pseudo-potential to the Vqo and V\\ terms for two iden- 
tical bosons and fermions, respectively, the eigenspectra for 
Vp S (r) and V m (r) agree very well for E < 0.5to (two iden- 



tical bosons) and E < 1.5 to (two identical fermions), and 
the unphysical eigenenergies are absent. This shows (i) that 
the scattering lengths aoo (identical bosons) and an (identical 
fermions) are dominant in this regime, and (ii) that the un- 
physical eigenenergies are due to the higher partial wave con- 
tributions of Vpg(r). The latter can be understood as follows: 
The coupling strengths gij'(k) for two interacting dipoles are 
proportional to ci[ji{k) /k ,+1 ' , and — since the a\ ji(k) are de- 
fined so that they approach a constant in the k — > limit — 
diverge as k goes to zero for / + /' > 0. A detailed analysis of 
the eigenequation for V vs (r) shows that these divergences give 
rise to the unphysical eigenenergies for small \k\. No unphysi- 
cal eigenenergies arise for larger \k\; in this regime, the 1 /k I+l 
factor in gu'{k) can be thought of as a simple "rescaling". 
Furthermore, the unphysical eigenenergies do not arise if the 
phase shifts are obtained for a short-range model potential 
whose scattering lengths are defined by — tan 8/ ji(k)/k l+l +1 . 
Although the pseudo-potential reproduces the eigenenergies 
well, we note that the single-parameter description fails to de- 
scribe the higher partial wave admixtures discussed in the con- 
text of Fig. |] 

In summary, this paper considers the scattering and bound 
state properties of two interacting dipoles near resonance. Al- 
though our analysis has been performed for a simple model 
potential, we believe that the main conclusions hold more gen- 
erally. Near resonance, the magnitude of all non-vanishing 
scattering lengths becomes large, with |aoo| being largest for 
identical bosons and an for two identical fermions. We have 
found that the wave function of weakly-bound two-dipole sys- 
tems contains higher partial wave contributions, raising in- 
teresting perspectives for studying BEC-BCS crossover-type 
physics. Despite the admixture of higher partial waves, a 
single-parameter pseudo-potential treatment reproduces the 
eigenenergy of the two-dipole system very accurately. 
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